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A note on the upward and downward intruder segregation in 
granular media 



Leonardo Trujillo 1 and Hans J. Herrmann 1 



Abstract The intruder segregation dependence on size 
and density is investigated in the framework of a hydro- 
dynamic theoretical model for vibrated granular media. 
We propose a segregation mechanism based on the differ- 
ence of densities between different regions of the granular 
system, which give origin to a buoyant force that acts on 
the intruder. From the analytic solution of the segregation 
velocity we can analyze the transition from the upward to 
downward intruder's movement. 



The understanding of the behavior of granular segrega- 
tion is relevant due to their practical importance in many 
industries fj],@|. When a granular mixture is subject to ver- 
tical vibrations under gravity, the grains tends to seg- 
regate with the larger particles at the top of the bed. 
Many studies have been devoted to this subject, usually 
referred as the "Brazil nut effect" |,|,g||,|,|,|,0,|ll|,0 



LJ,|FJ,|H||Hj,|l^,|lg. There exists a great controversy con- 



cerning the upward to downward segregation in granular 



materials |9|,g0[|21 
lem using a recent 



23LI24]. Here we address this prob- 
y proposed model for intruder size seg- 
regation in dry granular media[p5[. This model charac- 
terizes the rise velocity of a large intruder particle im- 
mersed in a medium of monodisperse fluidized small par- 
ticles. In Ref.[^5| we have proposed a segregation mech- 
anism based on the difference of densities between differ- 
ent regions of the system, which gives origin to a buoy- 
ant force that acts on the intruder. This force include 
an Archimedean buoyancy force due to the differences be- 
tween the intruder material density pi and the bed den- 
sity pf, f.4 = (pf — Pi)Vig. Here g is the gravity field and 
Vj = —jj-rf is the D-dimensional volume of an intruder 
with radius r 7 . The factor Q D = 2n D/2 / r(D /2) is the sur- 
face area of a D-dimensional unit sphere. Also, we include 
a thermal buoyancy force caused by density variations of 
the granular fluid which comes from differences in the local 
"granular temperature" AT g . The change in the granular 
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fluid density through the thermal expansion, produced by 
the difference of temperatures, is p' F = /9f(1 — aAT g ), 
where a is the thermal expansion coefficient. The ther- 
mal contribution to the buoyancy force is f-r = AppVig, 
where App = p'p — pF = —appAT g . The granular temper- 
ature T g is defined proportional to the mean kinetic energy 
associated to the velocity of each particle. The granular 
temperature difference AT g is due to the dissipative na- 
ture of the collisions between grains. This difference is due 
to the fact that the number of collisions on the intruder 
surface increases with the size, but the local density of 
dissipated energy diminishes. The region with intruder is 
hotter than the region without intruder. This lead to the 
thermal buoyancy force that contributes to the intruder's 
upward movement. The intruder also experiences a viscous 
drag force of the granular fluid. The drag force id [^6[|2^] is 
considered to be linear in the velocity of segregation u(t), 
and is like the Stokes' drag force id — —6Trprju(t), where 
p is the coefficient of viscosity. Therefore the equation of 
motion that governs the segregation process is 



Qd d du(t) f2 



D 



■rj Pr 



dt 



^-rf [ Pl -p F (l + aAT g )} g 
-6nprju(t). 



(1) 



We take the reference frame positive in the upward 
vertical direction. Arranging terms in Eq.(Q) we find the 
following differential equation 



du(t) 
dt 



= [p F (l + aATg) - PI ] 



and the solution is 



a 



(t) = t g 



^(1 
Pi 



a 



AT„ 



Pi f2 D pirj 



6 ^(t), (2) 



where the time-scale to is 
to 



&DPir?~ 



(3) 



(4) 



6irDp 

The drag force always acts opposite to the intruder ve- 
locity. So, the intruder's upward/downward movement is 
exclusively due to the buoyancy forces. For our theoretical 
calculations we define the settling velocity u s 



to9 
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(5) 



When u s > the resulting movement is upward. On the 
other hand, if u s < the resulting movement is down- 
ward. Our analysis reveals that if (1 + aAT g )pp / pi > 
1, the intruder ascends. The opposite occurs when (1 + 
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(a) 



(b) 



Fig. 1. Regions used to calculate the granular temperature, 
(a) Region around the intruder within a sphere of radius ro 
and (b) region without intruder. 



uAT g )pF I pi < 1. When u s = there is no upward nei- 
ther downward movement. 

The intruder's presence modifies the local temperature 
of the system due to the collision that happen at its sur- 
face. In Ref. p5[ we have proposed an analytic procedure 
to estimate the temperature difference among the granu- 
lar fluid. We can calculate within a sphere of radius ro the 
value of the temperature T± in the granular fluid in pres- 
ence of the intruder and compare it with the temperature 
T2 in the granular fluid without intruder. In both cases we 
calculate the granular temperatures (TijTa) at a distance 
r = ri from the center of the sphere of radius ro. See Fig. 1 
for a schematic picture of the regions used to calculate the 
granular temperature. In the model this temperature ratio 
for two dimensions is given by (see Ref. pa] for details) 



Tx = / / o (A F r o )[0 Ai3 /o(A F r J ) + K (X F r 7 )] 
T 2 \Io{XFri)[0 AB Io(X F ro) + K (\ F r Q )} 



(6) 



where Iq(x) and Kq(x) are modified Bessel functions, and 
A F / (A/r / )X 1 (A F r / ) + Xih(Xjrri)K (.X F rj) 



e = 



(7) 



A F / (A/r/)/i(A F r/) - Xih(Xiri)h(X F ri) 
For three dimensions, 

Ti _ ( ^o(A F ro)[6)A5»o(A F rj) + fc (A F rj)] 
T 2 \«o(A F r/)[6> J 4B«o(A F ro) + fc (A F r )] 

where io(x) and ko(x) are spherical modified Bessel func- 
tions, and 



(8) 



e = 



A F z (Ajrj)fci(A F rj) + A/Zi(Ajrj)fc (A F rj) 
A F i (A/r/)ii(A F r/) - A/ii(A/r/)i (A F r/) ' 



(9) 



The granular bed is formed of N monodisperse parti- 
cles of mass m F and radius r F . The particles are inelastic 
hard disks (D = 2) or spheres [D = 3). The inelasticity 

is specified by a restitution coefficient e < 1. The factor A 

1/2 

couples the coefficient of thermal conductivity k — KyTg 
and the dissipation rate 7 = £Tg^ 2 , explicitly one has 1 25 

s_ 

2k ' 



A 



(10) 



where £ and kq depend, among other things, on the mass 
of the particles. For the fluid particles this is 



C F = -^=(l-e 2 )n 2 ffo(2r F )^- 1 (J- 



1/2 



(11) 



and for the region near the intruder 

sD-l ( m i 
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,(12) 



where n = N/V, V = L D is the volume of the system 
of size L, go is the pair correlation function for two fluid 
particles. In 2D the pair correlation function is |2q| <?o = 
(1 — jqv) /(l — v) 2 1 with the area fraction v = n7rrjTln 3D 
the pair correlation function is[^9| go = (2 — v)/2(l — y) 3 , 
with the volume fraction v — Amrrp/3. Equations (|10[), 
( |TT|) and ( |l2| ) define the factors A F = ^/£ f /2k and A/ = 
yJ7/2ko. The prefactor kq in two dimensions is 

1/2 r 



Ko = 3nr F 
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and for three dimensions 
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,(14) 



where G = vgo. 

For a dense system the pressure is related to the den- 
sity by the virial equation of state is p — ^^nT g (l + 
Qg-ngo{2rp) D ). The thermal expansion coefficient is de- 
fined as a = y -1 (<9V '/ 'dT g ). From the equation of state we 
can calculate the coefficient a — T~ J (i/ 2 + 8) 2 / (v 3 — iv 2 — 
Sv — 8) (y — 1) for two dimensions and a — T^lv 3 — v 2 — 
v — l){v—l)/(v i — Av 3 + Av 2 + 4:V+l) for three dimensions. 

The state-dependent viscosity possesses the general 
form p — /io^/X^. The prefactor /xq is 



—nrp (amp) 1 1 2 



for two dimensions and 



— nrp^mp) 1 ! 2 



5 1 4 
1+ 16G + 5 



1 



— 1 o 

7T 



(15) 



(16) 



for three dimensions. 

The temperatures of the region with intruder and the 
region without intruder are different. For T1/T2 > 1 the 
thermal buoyancy force favors the upward movement. When 
T1/T2 < 1 the thermal buoyancy force favors the down- 
ward movement. The nonlinear form of Equations (||) and 
(||) doesn't allow us to calculate analytically the explicit 
dependence on the mass ratio m//m F and the size ra- 
tio cj> = rijrp. Let us examine numerically the tempera- 
ture ratio T1/T2 as a function of the mass ratio m//m F 
for different values of <f>. We set the number of particles 
N = 5 x 10 3 , the volume fraction v = 0.75, r = L/2 
and the restitution coefficient e = 0.95. We consider small 
particles with unitary mass mp = 1, and we vary the in- 
truder's mass like m/ = xm Fl which x ~ (mj << m F ) 
to x = 2. Figure 2 shows the temperature ratio depen- 
dence as function of the mass ratio for two and three di- 
mensions. 

From the analysis of Fig. 2 we can conclude that al- 
ways Ti > Ti- and the upward to downward movement is 
basically controlled by the Archimedean buoyancy force. 

Let us concentrate on the dynamic of the segregation 
process described by the Equation (|J). From the settling 
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Fig. 2. Granular temperature ratio T1/T2 as a function of 
the mass ratio mi/rriF for (a) two dimensions, and (b) three 
dimensions. 



velocity u s (Eq. (JH])) we can analyze the dependence on 
size and density ratio. An uniformly fluidized state can be 
realized when the granular system is subject to a verti- 
cal vibration with amplitude Aq and frequency ujq = 2nf. 
In the experiments the excitation is described by the di- 
mensionless acceleration r = Aou>o/g. The characteristic 
velocity of the system is uq — Aqujq. The system increases 
its energy as a result of external driving while its de- 
creases its energy by dissipation. In our theoretical model 
we do not consider a sinusoidal excitation. Analogies with 
shaken granular systems and the dependency of the gran- 
ular temperature on the amplitude of vibration has been 
studied in Ref.j50| for dense granular systems in 2D, in 
which the following expression relating the global granu- 
lar temperature to a symmetric vibration with maximum 
velocity u Q is: T g = 2^m F L(A Lo ) 2 /2Nr F {\ - e 2 ). For 
three dimensions we estimate the granular temperature as 
T g ~ rriF {Aqloq) 2 . In order to calculate u s we use the fol- 
lowing model parameters: mass particle density pf = 2.7 
gem" 3 , r F = 0.1 cm, e = 0.9, v = 0.75, N = 5 x 10 3 , 
g = 100 cms -2 , r = L/2, A = 2r F cm, w = 5.81y/g/Ao 
for two dimensions, and us® = 0.7-^/ g/Ao for three dimen- 
sions. Fig. 3 summarizes the results of our calculations 
for size ratio 4> from 1 to 10 in two dimensions and three 
dimensions. 

From these results we can observe that the segregation 
is rapid in systems in which small particles are more dense 
than the intruder. For fixed values of (f>, when the density 
ratio increases the rise velocity diminishes, this implies 
a diminution of the intruder rise time. Finally, when the 
intruder's density is bigger than the small particles density 
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Fig. 3. Dimensionless intruder settling velocity u s / (Aoloo) , as 
a function of density ratio pi/pF in (a) two dimensions, and 
(b) three dimensions. 
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Fig. 4. Phase diagram determined for (a) two dimensions, and 
(b) three dimensions. 



the intruder particle sinks u s < 0. For a fixed value of the 
density ratio the segregation rate increases with the size 
ratio. These result shown the condition for the crossover 
from the upward to the downward intruder's motion. 

Recently, the upward /downward transition has been 
studied in references [El, 3 1,32]. In Rcf. pjj Hong et ah, 



performed molecular dynamic simulations, in two and three 
dimensions, of weakly dissipative particles under gravity 
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that are in global thermal equilibrium with a heat reser- 
voir. They observed that for a binary mixture of granu- 
lar particles the large particles can rise and the upward/ 
downward movement depends on mass and diameter ra- 
tios. They proposed an explanation based on a compe- 
tition between the percolation effect and the condensa- 
tion of hard spheres under gravity (See Refs.|2l],|33|,[34| 
for a detailed explanation of the percolation-condensation 
mechanism) . Both and Hong presented a theory based on 
the variational principle for hard spheres and disks un- 
der gravityj3l]]. They also characterized the segregation 
phenomenon and investigated the crossover between the 
upward/downward movement for large particles and the 
dependence on mass and size ratios. In the framework of 
kinetic theory for a binary granular mixture, Jenkins and 
Yoon obtained a segregation criteria for spheres and disks 
that differ in size and/or mass[^2). Their mechanism is 
based on a competition between the inertia of the par- 
ticles through the ratio of partial pressures. The three 
different approaches presented in references [^l],|3l],|3^] co- 
incide qualitatively among them. In the same spirit we 
can calculate a phase diagram for the crossover from the 
upward/downward movement (See Figs. 4(a) and 4(b)). 
In our model these phase diagrams are derived from the 
dependence of the segregation velocity on the mass and 
size ratio at the situation where the intruder tends to rise 
(u s > 0), or to sink (u s < 0) (Figs. 3(a) and 3(b)). 

In summary, we have studied the segregation depen- 
dence on size and density for a single intruder in a fluidized 
bed. On a qualitative level our model is in agreement with 
the experimental phenomenology described in Refs.j|,[I]|: 
the higher the density of the intruder the lower the ten- 
dency to rise. In Fig. (3) we can note that when the density 
of the particles increases, in the range where pi/pf < 1, 
the segregation velocity diminishes, so the intruder's rise 
time will be bigger. In general a lighter particle rise more 
quickly than a heavier particle of the same size in the range 
where pi/pf < 1- In the range where pi / pf > 1 the in- 
truder sinks, and in this case the downward segregation ve- 
locity increases with the density and size. Our model also 
agrees with the qualitative behavior reported in Ref . p"7| , 
where the intruder height decreases as the density ratio 
increases. In this case, for the range where pi/pF < 1, the 
height of a lighter particle will be bigger than the height of 
a heavier particle of the same size at the same time. This 
is due to the fact that the rise velocity of a lighter particle 
will be bigger than the rise velocity of a heavier particle 
of the same size. In agreement with the theoretical results 
obtained in Refs.j2l],|[l],[32| and the numerical simulations 
of Refs. [p/7[|2l]| , our model predicts that in the case where 
the density ratio increases, the intruder particle will sink. 

A direct comparison with experiments performed in 
Ref. |25| ] are not appropriate in this work since they inves- 
tigate the intruder's density effect in presence of convec- 
tion and interstitial air. Comparison is also difficult with 
the experiments performed in Refs.[ ^9|[2^ |, where the role 
of inertia and interstitial air may play an important role. 
It is important to note that the experiment performed by 
Shinbrot and MuzzionJl9|], doesn't correspond with the 
typical experimental conditions of the Brazil nut effect. In 
these experiments the intruders are placed at the surface 



of a vibrated bed. They reported that the heavy intruder 
remained at the surface while the lighter intruder sank. 
The boundary condition for an intruder in this experi- 
mental set-up is more different and complicated than the 
boundary condition for an intruder immersed in the gran- 
ular bed. So, in this experiment we can not consider up- 
ward/downward transition. In the experiment performed 
by Liffman et ai, G2| they reported that the total rise 
time of the intruder is inversely proportional to the den- 
sity of the intruder. This is, a priori, in contrast to our 
findings. In our case, if the segregation velocity decreases 
with the intruder's density, then the rise time decreases. 
This apparent contradictory fact can be understood from 
the procedure followed by Liffman et ai, for the intruder's 
density variation. They studied the motion of disks placed 
at the base of a vibrated granular bed. They varied the 
density of the disks with similar size, drilling holes in the 
disks. The intruder is "transformed" from a disk to a ring. 
The inertia of these two object is different. Also, it is im- 
portant to note that the energy dissipation is different. 
Energy is dissipated during collision, among other things, 
due to the excitation of the internal modes of the ob- 
ject. So, the coefficient of restitution should be different. 
In terms of our model, based on the concept of granular 
temperature, the density of dissipated energy should be 
bigger for a ring than a disk. In this case the disk is "hot- 
ter" than a ring, and the buoyancy force would favor the 
upward movement. 

In order to compare our theory to the experiments one 
should do measurements varying the intruder's density ho- 
mogeneously in dry granular materials. Strangely, a sys- 
tematic study has not get been reported in the literature 
and since this should be easy to carry out, we hope that 
some experiment and/or simulation will soon be done. 
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